Electrostatic confinement of electrons in an integrable graphene quantum dot 
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We compare the conductance of an undoped graphene sheet with a small region subject to an 
electrostatic gate potential for the cases that the dynamics in the gated region is regular (disc- 
shaped region) and classically chaotic (stadium). For the disc, we find sharp resonances that narrow 
upon reducing the area fraction of the gated region. We relate this observation to the existence 
of confined electronic states. For the stadium, the conductance looses its dependence on the gate 
voltage upon reducing the area fraction of the gated region, which signals the lack of confinement 
of Dirac quasiparticles in a gated region with chaotic classical electron dynamics. 

PACS numbers: 73.63.-b, 73.63.Kv, 73.23.-b 
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A characteristic feature of the massless Dirac electrons 
in graphene, a single layer of graphite, is the suppression 
of backscattering [l| and the resulting strong angle de- 
pendence of scattering at interfaces [2| (related to "Klein 
tunneling" For clean graphene, graphene without 

atomic-scale defects or impurities that cause intervalley 
scattering, the transmission probability through an elec- 
trostatic barrier is unity at perpendicular incidence and 
drops sharply upon increasing the angle of incidence (cf. 
Fig. QJl) . The absence of backscattering at perpendicular 
incidence is a consequence of a combination of an effec- 
tive time reversal symmetry and band topology 
The same combination is responsible for the fact that, 
generically, Dirac electrons can not be confined in electro- 
statically defined structures or by a smooth random elec- 
trostatic potential through Anderson localization 0, 0| • 
For this reason, most experimental [H, 0] and theoreti- 
cal 

EES 

studies of confined electrons in graphene are 
performed on etched structures, in which inevitably in- 
tervalley scattering is strong and the physics essentially 
reduces to that of normal semiconductor quantum dots. 

This general argument against the possibility to con- 
fine electrons electrostatically overlooks a special case: 
Dirac electrons can be confined electrostatically if their 
classical dynamics inside the gated area is integrable and 
the corresponding Dirac equation is separable. The inte- 
grable classical dynamics allows for the existence of elec- 
tron paths for which the electron never approaches the 
boundary of the confinement area at perpendicular in- 
cidence. For such paths, an electrostatic barrier is fully 
reflecting, so that the electrons arc effectively confined. 
In quantum mechanics, it is the corresponding separabil- 
ity of the Dirac equation that invalidates the topological 
arguments against the electrostatic confinement of elec- 
trons @,[I3. 

The requirement of integrable dynamics is essential for 
the electrostatic confinement of the electrons. Certainly, 
if the classical dynamics is chaotic, electrons can not be 
confined electrostatically. There is a rich history of ex- 
amples in which the question of intcgrability of the clas- 




Figure 1: Angle-dependent scattering of an electrostatic inter- 
face (a), a circular quantum dot (defined by an electrostatic 
gate potential) in an infinite graphene sheet showing an ex- 
ample of a path that corresponds to a bound states (b), and a 
two terminal transport setup that can probe for the existence 
of bound states (c). 



sical dynamics has profound consequences for a quantum 
mechanical system. The best known example is that of 
the level spacing distribution of ballistic normal-metal 
quantum dots 1 



14j, which is described by Wigner- 



Dyson statistics for quantum dots with chaotic classical 
dynamics, and Poisson statistics for quantum dots with 
integrable classical dynamics. Another example, similar 
in spirit to the present study, is the lifetime of resonant 
states in optical resonators or microlasers [TBI, [l6| . Regu- 
lar lasers have optical ray trajectories that always exceed 
the critical angle and thus have an infinite lifetime, while 
chaotic trajectories are always eventually refracted out 
and thus have a finite lifetime fl7l ] . 

Experimentally, the existence of bound states would be 
clearly revealed in a two terminal conductance measure- 
ment setup schematically shown in Fig. [IF,. A "buffer" of 
undoped graphene separates a gated area from the leads 
and acts as an effective barrier. (The barrier is only "ef- 
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fective" because the conductivity of undoped graphene 
is 4e 2 /-7r/i, not zero.) A metal gate defines a "quantum 
dot" in which the carrier density is nonzero. If the dot 
is located approximately halfway between the metal con- 
tacts, eventual bound states in the dot will appear as res- 
onances in the conductance upon varying the gate voltage 
of the dot. Since only regular dots can support bound 
states, the presence of sharp resonances in the two ter- 
minal conductance gives a clear experimental signature 
of intcgrability in the graphene quantum dot. 

The remainder of this letter is reserved for material 
in support of the above scenario. We begin by consid- 
ering, as in Fig. [TJd, a circular quantum dot in an in- 
finite graphene sheet. This problem can be solved ex- 
actly, and the location and character of the bound states 
can be identified explicitly. We then turn to the two- 
terminal transport geometry of Fig. [T]: and compare nu- 
merical simulations of transport through a layer of un- 
doped graphene with circular and stadium-shaped quan- 
tum dots 

Taking the gate potential to be smooth on the scale of 
the lattice spacing, the electron dynamics in graphene is 
accurately described by the single valley Dirac Hamilto- 
nian 



H = v F p ■ a + V s&tc (x, y) 



(1) 



where a = {a x , a y ) are Pauli matrices and v-p is the Fermi 
velocity. For a circular dot we take the gate potential to 
be 

tWr) = hv F V d{r - R), (2) 

with $ the Heaviside step function. (The choice of a step 
function allows a closed-form solution; it is not essen- 
tial for the existence of bound states. The structure of 
quasi-bound states in the inverted setup [zero potential 
inside, nonzero outside] was considered in Ref. Il8l ) We 
look for eigenstates of Hq at zero energy. Because of the 
rotation symmetry of lg a t e , the eigenstates can be la- 
beled according to their angular momentum m, with m 
half-integral [l9(. They are 
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where the radial functions (p m ,tx with a — ± satisfy 

Vrn.a = -Vo&(R ~ r)<f m - a . (4) 
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For r > R the two equations (j4]) decouple and the solu- 
tions are 

¥W = a^- 1 ' 2 . (5) 

Normalizability requires \m\ > 3/2 and a+ = (a_ = 0) 
for m > (m < 0). For r < R one finds 

<Pm,+ (r) = bJ lm - 1/2 \(\Vo\r), (6a) 
<p m ,-(r) = -ibsgn(Vom)Ji m+1 /2\(\V \r). (6b) 



Continuity of the wave function at r = R then gives the 
implicit equation 



J\ m \-iM\V \R) = 0, 



(7) 



from which one obtains the values of the gate voltage Vo 
for which bound states exist at zero energy. These states 
correspond to states with nonzero orbital angular mo- 
mentum that circulate inside the disc, thereby avoiding 
Klein tunneling. The states are twofold degenerate, cor- 
responding to clockwise and anti-clockwise circulation. 
This twofold degeneracy is not specific to the circular 
shape of the quantum dot, but is rather an example of 
the Kramers degeneracy due to the effective time rever- 
sal, T = iijyC with C complex conjugation, of the single 
valley Dirac Hamiltonian 

The bound states exist for \m\ > 3/2 only. At angu- 
lar momentum |m| = 1/2 no normalizable bound state 
exists. Instead, the "zero mode" at \m\ = 1/2 is an ex- 
tended state. It is the existence of this extended zero 
mode that is responsible for the fact that electrons can 
not be confined electrostatically in a chaotic dot: For a 
chaotic dot, all wavefunctions will in general have a finite 
contribution from the zero mode, which will determine 
the finite lifetime of the state. This is to be contrasted 
with the case of the circle studied above, which admits 
solutions to the Dirac equation in which the zero mode 
has vanishing weight. 

Following Ref. 120 . transport properties in the two- 
terminal geometry of Fig. [TJ; are described by the Hamil- 
tonian 



H = H 



cad i 



(8) 



where Ui ea a(x) = for < x < L and Ui ea d(x) — > oo 
otherwise. The graphene sheet has width W. The con- 
ductance is obtained numerically using the transfer ma- 
trix method of Refs.[2ll andM. We place the center of the 
quantum dot at x = L/2. In Fig. [2] we show the results 
of such a calculation for the conductance for an undoped 
graphene sheet with a circular quantum dot, as a func- 
tion of the gate voltage Vq for aspect ratios W/L = 1 
and W/L — 6 and periodic boundary conditions in the 
transverse (y) direction. In both cases, we observe sharp 
resonances on top of a background conductance from the 
finite conductivity of the undoped graphene buffer. The 
positions of the sharp resonances are in good agreement 
with the gate voltages ([7]) of the bound states. We have 
also performed the same calculation with antiperiodic 
boundary conditions (results not shown) and obtained 
results that are qualitatively similar to those for periodic 
boundary conditions if W/L = 1 and indistinguishable 
from these if W/L — 6 (as was to be expected [20|). 

For aspect ratio W/L = 1, the background conduc- 
tance is essentially independent of the gate voltage Vq- 
For large aspect ratios, however, the background con- 
ductance has a relatively slow gate-voltage dependence 
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Figure 2: Conductance of an undoped graphene sheet with 
an electrostatically gated circular region of radius R (lower 
inset) placed in the center of the sheet, as a function of the 
gate voltage Vo for R/L — 0.26, W/L — 1 (upper panel) and 
R/L = 0.2, W/L — 6. The arrows show the positions of the 
solution to Eq. @ . Upper inset the dependence of the width 
of the first few resonances on dot size R/L. 



from a contribution from the zero mode. This is a fi- 
nite size effect. However, the density of states associated 
with the zero mode goes to zero only logarithmically in 
the limit L/R^> 1, which is why this feature persists for 
the entire range of system sizes accessible numerically. 
Still, the sharp resonances from the true bound states 
can be unambiguously identified, and their positions are 
independent of aspect ratio and boundary conditions. 

The width of the resonances in Fig. [2] depends on 
the angular momentum quantum number m of the res- 
onances. The dependence on m and on the size L of 
the undoped graphene buffer separating the dot from the 
electrodes can be understood in a simpler model that re- 
tains the circular symmetry. To that end we replace the 
leads by a circular electrode at a distance L from the 
center of the quantum dot. Since the integrated proba- 
bility of the bound states at distance L from the origin 
scales as (R/L) 2m , we conclude that the resonance width 
r should have the same dependence on L and m, 



r ~ (R/L) 



2\m\ 



(9) 



The larger the angular momentum quantum number m 
the narrower the resonances become. 



Figure 3: Same as Fig. [2] but with a stadium-shaped gated 
region (lower inset) with 2a/ R = y/H and R/L = 0.2f , W/L = 
f (upper panel) and R/L — 0.2, W/L = 6 (lower panel). The 
dependence of the sharpest resonance in the upper panel on 
the ratio R/L is shown in the upper inset. 



We have verified the prediction @ for the resonance 
width by repeating the numerical calculation for different 
values of R/L. The resonances have a Lorentzian shape, 
and we extract their width by fitting to the formula 
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(10) 



where xo is approximately given by the zero of the Bessel 
function (JT)). (There is a small shift of the resonance 
position from the coupling to the leads.) The dependence 
of the resonance width on R/L shown in the upper inset 
to Fig. [3] is in excellent agreement with the prediction ([9]) 
independent of boundary conditions and aspect ratio. 

Other regular shapes that give rise to integrable clas- 
sical dynamics and a separable Dirac equation behave in 
the same way as the circle. One analytically solvable ex- 
ample is given by a strip of width R extending across 
the entire width of the sample. For W ~ L the conduc- 
tance has separated Lorentzian resonances as a function 
of gate voltage, comparable to the circular gate voltage 
of Fig El reflecting discrete bound states. In the limit 
R <C L <c W these bounds states become effectively 
degenerate and their resonances start to overlap. The 
resonances are at position VqR = 7r(n + 1/2) with integer 
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Figure 4: Same as Fig. [2] but with an an electrostatically 
gated rectangular region of size R x 1.2R and rotated by 45 
degrees relative to the leads placed in the center of the sheet. 
The solid and dotted curves correspond to R/L — 0.2 and 
0.1, respectively. 



We can not numerically determine whether sharp reso- 
nances and, hence, confined electronic states, emerge at 
higher gate voltages, where the quantum-classical corre- 
spondence plays a more significant role. 

In conclusion, we have shown that gate potentials in 
which the (quantum) equation of motion is integrable, 
surrounded by undoped graphene, allow electrons to be 
confined. Gated regions in which the classical dynam- 
ics is chaotic, on the other hand, do not support bound 
states. Bound states manifest themselves in a two ter- 
minal conductance measurement through the presence of 
sharp resonances. Compared to the complexity of exper- 
iments on etched quantum dots in graphene, the setup 
suggested here is a promising candidate for the observa- 
tion of signatures of integrable dynamics in graphene. 

We are grateful to C. Kane for a stimulating discussion. 
This work was supported by the NSF under Grant No. 
DMR 0705476. 



n and their height diverges (in the limit U\ ea d — > oo), 
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cos VqR 



(11) 



As a prototypical example of a chaotic quantum dot we 
consider the conductance through a stadium (see lower 
inset in Fig. [3|). We have given arguments for why the 
stadium can not support bound states and we thus do not 
expect to see any sharp resonances in the conductance. 
This is confirmed by the calculated conductance trace in 
Fig. [3j which shows wide features which have a height 
that decreases and a width that remains constant upon 
increasing L/R (provided L/R is large enough). This is 
shown explicitly in the inset to the upper panel of Fig. [3] 
for the sharpest conductance peak near VqR — 6. Just as 
in the circular case contributions from quasibound states 
appear for larger aspect ratios. The behavior of these 
features is very similar to the zero mode contributions in 
the circle. 

For the Dirac equation, integrability of the classical 
equation of motion does not always imply integrability of 
the quantum mechanical equation of motion. An exam- 
ple is a uniformly gated region with a rectangular shape, 
for which the classical dynamics is integrable, but the 
Dirac equation is not [22j]. For the gate voltage range 
we considered in our numerical calculations (up to ~ 50 
levels in a circle of comparable size), we found no con- 
ductance resonances that have the same behavior upon 
increasing the ratio L/R as resonances in a circular quan- 
tum dot (becoming narrower at constant height); in con- 
trast all resonances in the range of gate potentials we 
considered become lower upon increasing L/R. (Typical 
data are shown in Fig.|U) This implies that, for this gate 
voltage range, a necessary condition for confinement is 
separability of the quantum mechanical equation of mo- 
tion, not merely integrability of the classical dynamics. 
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